De Haas-van Alphen effect in metals without inversion center 
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We show how the de Haas-van Alphen effect can be used to directly measure the magnitude of 
spin-orbit coupling in non-centrosymmetric metals, such as CePtaSi and LaPtaSi. 
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The recent discovery of superconductivity in a 
non-centrosymmetric heavy-fermion compound CePtsS: 1 
has renewed interest, both experimenta l 2 ^ 4 and 
theoreticalf^iSiL&^iiSiiiiiS to such materials. A peculiar 
property of non-centrosymmetric metals is that the spin- 
orbit coupling plays an essential role in the formation 
of single-electron states, namely it leads to the splitting 
of the energy bands characterized by helicity (i.e. the 
spin projection on the direction of momentum). This has 
important consequences for superconductivity: the elec- 
trons with opposite momenta have the same energies only 
if they are from the same non-degenerate band. For elec- 
trons from different bands this is possible only at some 
degeneracy lines or points in momentum space. There- 
fore, a large enough band splitting prevents the Cooper 
pairing of electrons from different bands. 

Theoretically, the magnitude of the band splitting can 
be determined from the band structure calculations. On 
the other hand, one can obtain some experimental in- 
formation about it from the frequencies of de Haas- 
van Alphen (dHvA) oscillations of magnetization. The 
first dHvA measurements in non-centrosymmetric met- 
als have been reported in Ref. While the restoration 
of the Fermi surface in CePt 3 Si is difficult due to large 
values of the effective masses, the measurements on its 
light-electron counterpart LaPtaSi have revealed rich in- 
formation about the band structure. 

Previous experimental work on the dHvA and a closely 
related Shubnikov-de Haas effects in systems without in- 
version center focused either on asymmetric semiconduc- 
tor heterostructureSfi^* 1 ^* 1 ^ or on bulk semiconductors 
with zinc-blende structured A common feature of these 
systems is that the spin-orbit band splitting results in two 
distinct frequencies of the dHvA oscillations. When the 
frequencies are close, their interference produces a char- 
acteristic beating pattern in the observed signal. This 
phenomenon was first theoretically predicted in Ref. 
(for recent work on the subject, see e.g. Ref. [Til) . An- 
alyzing the beating pattern allows one to estimate the 
strength of the spin-orbit coupling. In this brief article 
we apply these ideas to the interpretation of the dHvA 
data in CePtsSi and LaPt 3 Si. 

The effective single-electron Hamiltonian in a non- 
centrosymmetric crystal can be written in the form 

H = e (k) + a(k)cr - fi B Hcr, (1) 

where eo(fe) is the band energy, the spin-orbit coupling is 



described by a pseudovector function a(k) = —a(—k), 
and er = (a x ,a y ,a z ) is the vector composed of Pauli 
matrices. The last term describes the Zeeman interac- 
tion with an external magnetic field H, with fiB being 
the Bohr magneton [using a general form of the Zee- 
man energy for band electrons, Hij{k)Hi<jj, would not 
add anything to the substance of our results] . The or- 
bital effect of the field can be included by replacing 
k — > k + (e / he) A(f) r& where f = is the position op- 
erator in the fc-representation and e is the absolute value 
of the electron charge. 

The momentum dependence of the pseudovector a(k) 
is determined by the point symmetry of the crystal. In 
the case of the tetragonal group C4 V , which describes 
the symmetry of both CePtsSi and LaPtaSi, it can be 
written quite generally in the form a(k) = a±[ipE(k) x 
z] + a z tpA 2 (k)z, where ifE and ipA 2 transform according 
to the irreducible representations E and A2 respectively, 
and a± and a z are constants £ The simplest polynomial 
expression compatible with the symmetry requirements 
is 



a(k) = a^(k y x - k x y) 



rk x kyk z (k x 



kl)z. (2) 



Setting a z = here we recover the Rashba model, 21 
which is used to describe the effects of the absence of 
mirror symmetry in semiconductor quantum wells. In 
cubic zinc-blende crystals, the momentum dependence of 
oc(k) is given by the so-called A; 3 , or the Dresselhaus, 
term. 22 : 23 

One cannot expect Eq. Q to fully reproduce the spin- 
orbit band splitting in CePtsSi and LaPtsSi, which have 
quite complicated, multi-sheet, Fermi surfaces. Never- 
theless, this expression already captures the most im- 
portant, symmetry- related, features of the spin-orbit 
coupling, including the qualitative difference in the fc- 
dependences of a x ^ y {k) and a z (k), the presence of a band 
degeneracy line at k x = k y = 0, and the vanishing of 
a z (k) in the high-symmetry planes. A natural question 
is whether one can determine the strengths of both the 
xy- and z-components of the spin-orbit coupling using 
the dHvA experiments. 

The eigenvalues of the Hamiltonian JI} are 



e x (k) = e Q (fe) + A|a(fc) - ^ B H\ 



(3) 



where A = ± is the band index (note that the energy 
bands are split even at H = 0, if the spin-orbit coupling 
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is non-zero). There are two Fermi surfaces determined 
by the equations 



e\(k) = e F , 



(4) 



where e_F is the Fermi energy. Although there may be ac- 
cidental degeneracies at some magnitudes and directions 
of the field, in general there are no symmetry reasons for 
the Fermi surfaces to intersect. Indeed, this would hap- 
pen if a(k) = HbH. These three equations can have so- 
lutions at some isolated points in the first Brillouin zone, 
which may or may not be on the Fermi surface. The 
shape of the Fermi surfaces (QJ depends on the magnetic 
field, which can be directly probed by dHvA experiments. 
In particular, while at H = we have e x (— k) — e\(k), 
which is a consequence of time reversal symmetry, in the 
presence of magnetic field the time-reversal symmetry is 
lost, and e x (— k) ^ e\(k), i.e. the Fermi surfaces do not 
have inversion symmetry, in general. 

To calculate the dHvA frequencies, one needs to in- 
clude the coupling of the magnetic field to the orbital 
motion of electrons. In the quasi-classical approximation 
one can derive the Lifshitz-Onsager quantization rulesj^S 
which implicitly determine the energy levels of the band 
electrons: 



S\(e, k H ) 



2ireH 

he 



[n + 7A (r)]. 



(5) 



Here S\ is the area of the quasi-classical orbit, T, in the k- 
space defined by the intersection of the constant-energy 
surface e\(k) = e with the plane k ■ H = kn (H 
11 11 ) . n is a large integer number, and < 7a (r) < 1 is 
a constant, which depends on the Berry phase acquired 
by a band electron as it moves along FiSiSi The value 
of 7a (r) does not affect the expressions for the dHvA 
frequencies discussed below. 

The oscillating magnetization contains contributions 
from both bands and can be approximately written as 



M n 



A 



M\ cos 



( 2ttF x 
\ H 



(0) 



where M x and <p x are the amplitudes and phases of the 
oscillations. The expressions for the amplitudes are given 
by the standard Lifshits-Kosevich formulas i^S The dHvA 
frequencies F x are related to the extremal, with respect 
to fcff , cross-sectional areas of the two Fermi surfaces as 
follows 



he 
2ne' 



St 



(7) 



[in addition to the fundamental harmonics I©, the ob- 
served dHvA signal also contains higher harmonics with 
frequencies given by multiple integers of Fa]. 

If the external field is weak compared to the spin-orbit 
band splitting, i.e. \±bH -C |«(fe)|, the band energies © 
can be represented as a Taylor expansion 



e\(k) = eo(fe) + A|a(fe)| - A^ s (di2~) 
[H 2 -(aH) 2 } + ..., 



A/z 2 



B 

~%a.ik\ 



(8) 



where a(k) — a(k) /\a(k)\. Similarly, the extremal 
cross-section areas can be written in the form 



S e x xt (H) = S e x x \0) + A X {H)H + B X {H)H 2 



(9) 



The second, linear in H, term on the right-hand side pro- 
duces the phase shifts in dHvA signal This effect is 
similar to the usual phase shift due to a paramagnetic 
splitting of Fermi surfaces in ccntrosymmctric metals. 
For some directions of the field, the linear term can be 
absent, see an example below. The third term and all the 
subsequent terms produce the magnetic field dependence 
of the dHvA frequencies. This is a specific feature of the 
dHvA oscillations in crystals without inversion symme- 
try, which can be observable if the Zeeman energy is not 
too small in comparison with spin-orbit coupling. A non- 
linear field dependence of the dHvA frequencies has been 
observed in asymmetric quantum wells 

To illustrate the above statements, let us look at a 
simple example of a three-dimensional elliptic Fermi sur- 
face with eo(k) = h 2 k 2 _/2m± + h 2 k z /2m z — ep, where 
k±_ = (k x , k y ), and m±, m z are the effective masses. The 
Fermi momentum kp is introduced via ep = h 2 k 2 F /2m±. 
We consider only H || z to make connection with the 
experimental results of Ref. 1131 where two main dHvA 
branches, named a and /3, were detected for this field 
orientation. One can show that the linear in H terms in 
the expansions JSJ and 10 vanish. The maximum cross- 
sections of the Fermi surfaces correspond to k z = 0, then 
l PA 2 = and we obtain the extremal cross-section area 
which depends only on the transverse spin-orbit coupling: 



S e x xt (H) 



Tiki 



,,2 172 

2a\k 2 F 



(10) 



In obtaining this result we used the expression J2J) for a 
and assumed that the Zeeman energy is small compared 
to the spin-orbit band splitting, which in turn is much 
smaller than the Fermi energy: fisH ^ |a^|/ci? <C £f- 
Although, for a more complicated Fermi surface, there 
might be additional extremal cross-sections at nonzero 
k z r^ the linear in H term in Eq. is still absent due 
to the symmetry properties of ifA 2 ■ 

To estimate the magnitude of the effects under consid- 
eration, we use the expressions Q and (fTUf> to calculate 
the difference of the dHvA frequencies: 



F_ - F, 



2c 
he 



a±\kpin± 1 



2a\k 2 F 



(11) 



The experimental measurement of the splitting of the 
frequencies allows one to determine the strength of the 
spin-orbit coupling. Using as an example the frequencies 
of the a and (3 branches from Ref. Eli 7 ^ = 1.10 x 10 8 Oe 
and Fp — 8.41 x 10 7 Oe, and mj_ ~ 1.5m, we ob- 
tain for the spin-orbit splitting of the Fermi surfaces: 
lajjfcf — 10 3 K. While the results of the band struc- 
ture calculations for LaPtaSi reported in Ref. [RJ do not 
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contain explicit values of the band splitting AE so , for 
CePt 3 Si one has AE so ~ 50 - 200meV£ As for the mag- 
nitude of the magnetic field dependence of the frequency 
splitting, in the range of fields used in Ref. [l3| (up to 
17T), we have ^ B H/\a^\k F ~ 1CT 2 . 

We would like to note that the expansions JHJ) and 
therefore 10 fail if ct(fe) = 0. According to Eqs. J2J, 
this happens if the extremal orbit passes through the 
poles of the Fermi surface, where the bands are degen- 
erate. In this case, the so-called "magnetic breakdown" 
occurs, in which the electrons can tunnel from one band 
to another near the degeneracy points. Instead of Eq. 
(|q]) , the dHvA signal then contains additional fundamen- 
tal harmonics corresponding to the quasi-classical orbits 
switching between different bandsi^l It is not clear if this 
phenomenon occurs in LaPtaSi and CePtsSi. 

In conclusion, we have discussed how the absence of 
inversion symmetry in the crystal lattice of a metal man- 
ifests itself in the dHvA experiments. The splitting of 
the dHvA frequencies is a direct measure of the parame- 
ters of the effective spin-orbit Hamiltonian. In particular, 



according to Eq. , it allows one to estimate the mag- 
nitude of the "transverse" component of the spin-orbit 
coupling (in contrast, there seems to be no simple way to 
determine the z-axis component using the dHvA data). 
Also, the interplay of the Zeeman and the spin-orbit in- 
teractions results in a deformation of the Fermi surface, 
which is responsible for a non-linear field dependence of 
the dHvA frequencies, the effect absent in centrosymmet- 
ric metals. 
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